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Abstract 


Considering  the  exploitation  needs  associated  with  the  Synthetic  Aperture  Radar  (SAR) 
applications  involving  moving  and  non -stationary  targets,  a  fundamental  spectral  domain  model 
for  moving  point  and  distribution  of  scatterers  is  presented.  The  approach  is  accurate  as  the 
spherical  phase  front  is  rigorously  treated  throughout  the  Doppler  analysis.  Due  to  the  analytic 
nature  of  the  model  in  range-Doppler  plain,  large  squint- angle  operations  (e.g.,  SpotSAR)  can  be 
characterized.  Furthermore,  motion  characteristics  can  be  extracted  in  the  sub- aperture  level  to 
better  reflect  the  general  motion  (e.g.,  non-uniform).  As  a  result,  enhanced  adaptivity  can  be 
achieved  to  improve  the  imagery  and  other  target-signature  measured  products.  The  range- 
Doppler  model  can  serve  as  a  tool  for  polari metric  and  multi-channel  analysis. 


Resume 


Compte  tenu  des  besoins  d' exploitation  associes  aux  applications  SAR  (radar  a  synthese 
d’ouverture)  utilisees  avec  des  cibles  mobiles  et  non  stationnaires,  nous  presen  tons  ici  un  modele 
fondamental  dans  le  domaine  spectral,  applicable  aux  points  mobiles  et  a  la  distribution  de 
diffuseurs.  Cette  approche  est  exacte,  car  le  front  de  phase  spherique  est  traite  rigoureusement 
tout  au  long  de  T  analyse  Doppler.  Etant  donne  la  nature  analytique  du  modele  dans  le  plan 
distance-Doppler,  de  grandes  operations  a  angle  de  strabisme  (p.  ex.  SpotSAR)  peuvent  etre 
caracterisees.  En  outre,  les  caracteristiques  de  mouvement  peuvent  etre  extraites  au  niveau  sous- 
ouverture  pour  rnieux  representer  le  mouvement  general  (p.  ex.  non  uniforme).  Par  consequent, 
l'adaptivite  peut  etre  accrue  pour  ameliorer  Timagerie  et  d'autres  produits  mesures  de  signature 
de  cibles.  Le  modele  distance-Doppler  peut  servir  d’outil  pour  l'analyse  polarimetrique  et  multi¬ 
canal. 


DRDC  Ottawa  TM  2007-289 


This  page  intentionally  left  blank. 


DRDC  Ottawa  TM  2007-289 


Executive  summary 


A  basic  Fourier  transform  pair  for  slant  range-Doppler  modeling 
of  moving  scatterers  for  SAR  applications:  theory 

Ramin  Sabry;  DRDC  Ottawa  TM  2007-289;  Defence  R&D  Canada  -  Ottawa; 
November  2007. 


Background:  Synthetic  aperture  radar  (SAR)  provides  means  for  fine-resolution  detection  and 
imagery  by  virtue  of  the  synthetically  created  coherent  gain.  Conventional  SAR  technologies  use 
certain  established  techniques  to  characterize  and  exploit  the  coherent  phase  and  create  imagery. 
All  mainstream  techniques,  however,  use  phase  simplification  strategies  (e.g.,  plane- wave 
approximation  in  a  tomographic  paradigm)  to  approximate  the  spherical  phase.  Such  widely  used 
approximations  become  less  adequate  and  more  problematic  when  exploiting  the  moving  targets 
imagery  for  motion  modeling,  modifications  (e.g.,  focusing)  and  target  recognition.  An  analytic 
range-Doppler  model  to  characterize  canonical  moving  point  scatterers  is  developed  here.  The 
approach  is  accurate  as  the  spherical  phase  front  variations  and  Doppler  effects  due  to  the  target 
motion  are  directly  modeled  in  slant  range-Doppler  domain  without  the  need  for  range  and 
azimuth  de-coupling. 


Results:  The  system  response  or  Green’s  function  for  imaging  a  collection  or  distribution  of 
moving  scatterers  is  derived  and  established  by  taking  advantage  of  the  developed  coherent 
range-Doppler  model.  Application  of  the  range-Doppler  system  model  in  conjunction  with  the 
measured  data  and/or  image  yields  target  velocity  information.  On  this  basis,  a  velocity  extraction 
algorithm  is  also  developed  to  characterize  the  motion  of  the  scatterer  or  target.  The  range- 
Doppler  model  can  also  be  used  to  define  the  motion  corrected  matched-filter  transfer  function 
for  re-focusing  the  SAR  image. 


Significance:  Certain  classes  of  applications  would  largely  benefit  from  the  accuracy  of  the 
developed  technique  where  methods  based  on  the  conventional  phase  expansions  (e.g.,  series, 
Fresnel)  would  not  be  effective.  Squint  mode  of  operation  (larger  squint  angles,  in  particular)  can 
be  modeled  more  correctly  using  the  developed  approach.  Also,  wider  synthetic  aperture 
operations,  i.e.,  large  integration  time,  can  be  handled.  The  coherent  nature  of  the  model  offers 
accurate  motion  extraction  through  effective  utilization  of  the  phase  information.  Means  for  sub- 
aperture  analysis  is  provided  which  enables  characterization  of  the  motion  variations  within  the 
synthetic  aperture  acquisition  time. 


DRDC  Ottawa  TM  2007-289 


Sommaire 


A  basic  Fourier  transform  pair  for  slant  range-Doppler  modeling 
of  moving  scatterers  for  SAR  applications:  theory 

Ramin  Sabry;  DRDC  Ottawa  TM  2007-289;  R  &  D  pour  la  defense  Canada  - 
Ottawa;  Novembre  2007. 


Introduction  :  Le  radar  a  synthese  d'ouverture  (RSO)  permet  la  detection  et  l'imagerie  a  haute 
resolution  au  moyen  d’un  gain  coherent  cree  par  synthese.  Les  technologies  RSO  classiques 
utilisent  certaines  techniques  etablies  pour  caracteriser  et  exploiter  la  phase  coherente  et  produire 
des  images.  Toutes  les  techniques  courantes  font  toutefois  appel  a  des  strategies  de  simplification 
de  phase  (p.  ex.  approximation  d'onde  plane  dans  un  paradigme  tomographique)  pour  s’approcher 
de  la  phase  spherique.  Ces  approximations  largement  utilisees  sont  moins  satisfaisantes  et  plus 
problematiques  lorsqu'on  exploite  l'imagerie  de  cibles  mobiles  pour  la  modelisation  du 
mouvement,  des  modifications  (p.  ex.  mise  au  point)  et  la  reconnaissance  de  cibles.  Nous 
elaborons  ici  un  modele  Doppler  analytique  en  fonction  de  la  distance  pour  caracteriser  des 
diffuseurs  ponctuels  mobiles  canoniques.  Cette  approche  est  precise,  car  les  variations  du  front  de 
phase  spherique  et  les  effets  Doppler  dus  au  mouvement  des  cibles  sont  directement  modelises 
dans  le  domaine  Doppler  a  distance  oblique  sans  necessiter  de  decouplage  en  distance  et  en 
azimut. 


Resultats  :  Le  modele  Doppler  elabore  a  distance  coherente  permet  de  determiner  la  reponse  du 
systeme  ou  la  fonction  de  Green  pour  l'imagerie  du  regroupement  ou  de  la  distribution  de 
diffuseurs  mobiles.  L’ application  du  modele  du  systeme  Doppler  en  fonction  de  la  distance  et 
l'utilisation  des  images  et/ou  donnees  mesurees  fournissent  de  l'information  sur  la  vitesse  des 
cibles.  De  meme,  un  algorithme  d' extraction  de  la  vitesse  est  etabli  pour  caracteriser  le 
mouvement  d'un  diffuseur  ou  d'une  cible.  Le  modele  Doppler  en  fonction  de  la  distance  peut 
aussi  servir  a  definir  la  fonction  de  transfert  a  filtrage  adapte  et  correction  selon  le  mouvement 
pour  refaire  la  mise  au  point  d'une  image  de  RSO. 


Portee  :  Certaines  classes  d’ applications  profiteraient  largement  de  la  precision  de  la  technique 
elaboree  lorsque  les  methodes  fondees  sur  des  expansions  de  phase  ordinaires  (p.  ex.  en  serie  et 
de  Fresnel)  s’avereraient  inefficaces.  La  methode  etablie  permet  de  modeliser  plus  correctement 
le  mode  de  fonctionnement  a  strabisme  (plus  grands  angles  de  strabisme,  en  particular).  Elle  se 
prete  aussi  au  traitement  des  ouvertures  synthetiques  plus  larges,  c.-a-d.  aux  longues  durees 
d’ integration.  La  nature  coherente  du  modele  procure  une  extraction  de  mouvement  precise  par 
l'utilisation  efficace  de  l’information  de  phase.  Des  moyens  d'analyse  de  sous-ouverture  peuvent 
caracteriser  les  variations  de  mouvement  durant  1' acquisition  d’ouverture  synthetique. 


IV 
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1  Introduction 


Synthetic  aperture  radar  (SAR)  provides  a  means  for  fine-resolution  detection  and  imagery  by  the 
virtue  of  a  synthetically  created  long  antenna  as  the  result  of  antenna  or  sensor  motion  [  1  ]  -  [2] . 
Different  modes  of  SAR  operation,  e.g.,  single  channel,  polarimetric,  enjoy  growing  interest  due 
to  unique  properties  of  a  SAR  based  operation  for  remote  sensing.  In  addition  to  the  high 
resolution  property  and  the  ability  to  operate  in  non-cooperative  conditions,  advantageous  SAR 
characteristics  can  be  best  represented  by  the  three  conservation  principles  [1];  conservation  of 
energy,  confusion  and  coordinates.  The  first  two  principles  highlight  the  radiometric 
characteristics  where  the  latter  represents  the  spatial  performance  of  the  SAR  system.  These 
attractive  SAR  properties,  however,  are  accompanied  by  their  own  sets  of  challenges. 

Considering  the  coherent  nature  of  SAR  operation,  proper  treatment  of  phase  is  highly  important. 
Judicious  phase  modeling  and  characterization  becomes  even  more  critical  when  the  scattering 
centers  are  not  stationary.  In  such  scenarios,  care  should  be  taken  in  approximation  or 
decomposition  of  the  spherical  phase  front. 


Current  SAR  technologies  (e.g.,  airborne  and  space  borne  SAR,  PolSAR,  1SAR)  use  established 
techniques  to  characterize  and  exploit  the  coherent  phase  and  create  imagery.  In  range -Doppler 
processing,  the  range  dimension  of  the  image  is  established  through  the  resolution  of  fast  time 
delays.  The  azimuth  dimension  is  established  by  the  Doppler  gradient  of  the  radar  signal  phase 
function  across  the  scene.  All  mainstream  techniques,  however,  use  phase  simplification 
strategies  (e.g.,  plane-wave  approximation  in  a  tomographic  paradigm)  to  approximate  the 
spherical  phase.  These  approximations,  which  can  yield  good  results,  depending  on  the  applied 
technique  and  imaging  application,  are  needed  to  develop  a  practical  imaging  system.  However, 
the  same  argument  is  not  valid  for  the  advanced  exploitation  of  SAR  imagery.  Such  widely  used 
approximations  become  less  adequate  and  more  problematic  when  exploiting  the  moving  targets 
imagery  for  motion  modeling,  modifications  (e.g.,  focusing)  and  target  recognition. 


Despite  the  existence  of  significant  theoretical  work  and  applied  techniques  in  the  literature  [3]- 

[9] ,  it  is  not  easy  to  find  a  robust  algorithm  or  methodology  to  resolve  the  exploitation  and 
focusing  difficulties  associated  with  real  moving  targets  (e.g.,  non-uniform  motion,  moving  multi  - 
scatterers,  ships).  This  is,  undoubtedly,  a  great  challenge  due  to  many  contributing  factors  and  is 
not  intended  to  be  entirely  resolved  here.  Rather,  the  intent  is  to  address  the  challenge  by 
investigating  the  impact  of  accurate  phase  modeling.  The  present  approach  is  based  on  the 
analytic  modeling  of  point  and  distributed  moving  targets  (i.e..  Green’s  function  formalism)  in  the 
range-Doppler  domain  through  aperture  and  sub -aperture.  The  principle  of  stationary  phase  (PSP) 

[10]  is  employed  that  is  accurate  for  typical  large  time-bandwidth  SAR  signals.  As  such,  the  need 
for  phase  decoupling  (range -azimuth)  is  removed.  In  fact,  the  majority  of  advanced  techniques 
(e.g.,  [6]-[9])  aim  to  account  for  the  range  curvature  and  compensate  or  correct  for  range 
migration,  residuals  and  high  order  effects  resulting  from  the  spherical  phase  decomposition.  On 
this  basis,  the  approach  here  can  provide  better  accuracy  for  high  squint  angle  and  long  aperture 
modes  of  operation.  Moreover,  due  to  Fourier  transform  and  PSP  properties,  the  conventional 
plane  wave  approximation  is  avoided.  This  shift- varying  property  improves  the  range-Doppler 
accuracy  of  the  moving  scattering  centers. 
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Range-Doppler  frequency  domain  analysis  is  suitable  for  multi-channel  and  polarimetric 
investigations  as  the  need  for  time  domain  interpolation  is  removed.  This  interpolation  to  co¬ 
register  non-simultaneously  acquired  channels  (e.g.,  cross  polarimetric  channels)  is  used  for 
coherent  target  classification  analysis  and  often  is  the  source  of  error  resulting  in  incorrect 
interpretations  and  exploitations.  For  instance,  a  stationary  target  may  be  assumed  moving  due  to 
de-correlation  of  its  polarimetric  scattering  vectors. 


An  analytical  framework  and  methodology  is  established  here  with  pointers  to  the  type  of 
applications  that  benefit  from  this  approach.  At  the  present,  the  experimental  scope  is  limited  to 
comparing  the  simulated  canonical  moving  scatterer  range-Doppler  characteristics.  This  core 
comparison  would  indicate  the  relative  precision  achieved  for  the  ultimate  product  exploitation 
and  focusing. 


2 
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2  Theory 


Assume  a  moving  scatterer  with  large  time -bandwidth  radar  return  signal  product.  Figure  1 
depicts  the  associated  SAR  configuration  and  moving  scatterer  or  target  geometry.  The  SAR 
sensor  and  moving  target  positions  are  denoted  by  P  and  P  ,  respectively.  One  can 
mathematically  represent  the  echoed  SAR  signal  by: 


S(x,  r )  =  A(x)  P(r  -  R(x))  exp(-j2k  R(x))  exp (~ja(r  -  R(x))2 ) 


Where, 

A(x) :  Azimuth  antenna  pattern, 

P(r) :  Range  pulse  envelope, 

a  :  Pulse  fm  rate  determining  the  bandwidth  and  resolution, 

2tt 

k  =  — :  Propagation  wavenumber,  A  being  the  operating  wavenumber, 

A 

and  x,  r  denote  azimuth  and  slant  range  directions  respectively.  One  should  note  that  x,  r  are 
not  mathematically  orthogonal  and  should  be  treated  so. 


Figure  1.  Moving  Target  SAR  Imaging  Geometry 


DRDC  Ottawa  TM  2007-289 


3 


Represent  the  range  history  function  for  a  scatterer  undertaking  uniform  motion  with  a  generic 
form: 


R(x)=  yjR02  +  m2(x-x0)2 


(2) 


where  R()  is  the  equivalent  slant  range  of  the  closest  approach  (orthogonal  to  azimuth),  the  range- 
Doppler  slope  m  is  a  dimensionless  factor  representing  the  normalized  (by  stationary  target  rate) 
Doppler  rate,  and  x(j  represents  the  equivalent  or  shifted  scatterer  azimuth  position  (Doppler 

center).  For  a  stationary  target  (m  =  1 ),  R0  and  x0  become  the  conventional  slant  range  of  the 
closest  approach  and  scattering  center  azimuth  position.  The  range  history  function  in  (2)  can 
model  any  type  of  target  motion.  As  will  be  addressed  later  in  this  section  and  analyzed  in  the 
next  section,  general  motion  is  characterized  by  introducing  time- dependent  m  ,  R0  and  x0 ,  i.e., 

m(x) ,  R0(x)  and  x0(x) ,  in  the  range  history  equation  (2).  For  the  range- Doppler  variations  in 
this  section,  however,  a  constant  velocity  motion  is  assumed,  i.e.,  time -independent  m  , 

Rq  and  xQ  .  Employing  the  principle  of  stationary  phase  (PSP),  the  range- azimuth  frequency  (i.e., 
Doppler)  transform  of  the  waveform  (1)  is  calculated.  Due  to  non-orthogonality  of  the  range  and 
azimuth  pairs  (i.e.,  range  curvature),  conventional  direct  azimuth  transformation  neglecting  the 
curvature  is  not  analytically  correct.  Such  approximations  limit  and/or  suppress  insights  into  the 
exploitation/imaging  process  and  may  result  in  considerable  error.  To  resolve  this  difficulty  and 
maintain  the  analytic  consistency,  a  two  step  approach  is  undertaken  [5].  As  the  first  step  (that 
can  be  viewed  as  a  two-step  process  itself),  a  two  dimensional  correlated  Fourier  transform  is 
performed.  This  is  done  by  performing  the  range  Fourier  transformation  followed  by  an  azimuth 
Fourier  transformation  when  the  azimuth  dependence  of  the  range  curvature  is  taken  into 
considerations.  The  second  step  would  be  to  perform  an  inverse  range  Fourier  transformation  to 
the  range -Doppler  frequency  domain.  The  described  process  can  be  summarized  by  the  schematic 
depicted  in  Figure  2. 


To  apply  PSP  for  Fourier  transformation,  one  must  solve  for  the  stationary  point  of  the  total  phase 
associated  with  the  Fourier  integrand,  i.e., 


^  total  =0 


(3) 


and  use  the  result  to  compute  the  Fourier  integral  in  the  vicinity  of  the  stationary  point. 

dxJfX’r(x,  r) 

Employing  PSP  in  the  slant  range  ( r )  dimension,  i.e.,  - ’ —  =  0  ,  one  gets  (according  to 

dr 


(1)): 


4 
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S(x,ujr)  =  C0A(x)P(-^)exp(jU^)exp(-j(2k  +  ujr)R(x))  (4) 

2  a  4a 


where  'jJr  (rad/m)  is  the  range  space  frequency. 


Re-employing  PSP  in  the  azimuth  (x)  dimension,  i.e., 


dVx'u'(x,vr): 

dx 


stationary  point  equation  (according  to  (4)): 


0  ,  one  must  solve  the 


d  R(x)  'jJx 

dx  x=x°--(2  k  +  cor) 


(5) 


where  'jJx  (rad/m)  is  the  azimuth  space  frequency.  Using  (2)  one  obtains: 


%,ur)  =  C,%- 


t)- 


2km2[(l  +  ^f 

2k  4k  m 


P(-^)exp(j^\. 
2  a  4a 


(6) 


exp(-j2kR0[(l+^-)2-  W 


2k  4  k2~2 


]2)exp  (-ju}xxb) 


m 


In  order  to  derive  the  range-Doppler  characteristics  of  the  SAR  signal,  an  inverse  range  Fourier 
transform  is  applied.  This  would  require  solution  of  the  stationary  point 


equation 


&*“**>' (ux,ur) 
du)r 


0,  which  takes  the  form  (according  to  (6)): 


2 a 


^a+^> 


2  I 

U>r  ^2  ^ x  T9 

'  .  -  O  O  J 


(7) 


[(1 


2k  4k2  m2 


Expanding  the  term  associated  with  the  square  root  denominator  in  (7)  using: 
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1 


(8) 


(1  +  6)  2  =  1 


1  3  9  9 

— 6  +  -62+0(63) 

2  8 


and  assuming  a  typical  bandwidth  (linewidth): 


'■°r  ,  ^r2 

k  4  k2 


1 


CJ, 


a,  2  2 

4  k  m 


<1 


(9) 


one  can  rewrite  (7),  for  the  phase  to  the  order  of  0(cj  ,3 )  as: 


^y+r-Rn(u;x)=  0  (10) 

2a 


where 


a'  a 


*0  ,jJx 


4k  3m  2(  I  - 


4fc  2ffZ 2 


(1 


_*0 _ 

2  1 
_ 9  2 

4k2m2 


(ID 


(12) 


Applying  the  solution  of  (10)  for  the  stationary  point  jj  r  and  using  (6),  one  derives  the  range - 
Doppler  description  of  the  SAR  signal  as: 


6 


DRDC  Ottawa  TM  2007-289 


SK.r)  =  C,  A(x0 - - - - — )• 

lkm2[{l- S-fr-i^K)))2  -  -%—P 
k  4k  m 

P(  —  (r-Rn(ux)))exp(j  —  (r-Rn(ux))2)-  (13) 

a  a 

/  2  I 

exp(-j2kR{}  [(1-  (r~  Rn (ojx)))2  —  ^  \2  )exp(-  juJxx0) 

k  4k  m 


C0 ,  C, ,  C2  in  (4), (6)  and  (13)  are  complex  constants. 

The  above  equations  suggest  that  one  could  ultimately  define  a  range-Doppler  domain  Green’s 
system  function  for  multi -scatterers  as: 

G(  r,  u>x,  X,  Y,  m)  =  S(ux,  r)  (14) 

where  X ,  Y  determine  the  initial  position  of  the  moving  scatterer.  One  can  express  the  target 
equations  of  general  motion  as: 


X(t)  =  x  +  Vx(t)t 

(15) 

Y(t)=Y+Vy(t)t 

(16) 

where  Vx(t)  and  Vy(t)  are  the  generalized  velocity  functions.  Considering  the  sensor  or  radar 

x 

platform  altitude  H  and  velocity  va  ,  and  setting  t  —  —  one  can  write  the  general  range  history 

vfl 

function  of  a  moving  scatterer  as: 

R2(x)  =  (x  -  X (t))2  +  Y (t)2  +  H2 

=  (x- X -7r(x)x)2  +  (F  +  jy(x)x)2  +  H2 


where 
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(18) 


V, 


7*00  = 


7vOO 


V 

X 

O7 

,v 

(19) 


Expanding  (17)  and  comparing  with  (2),  for  general  motion  one  obtains: 

m(  x)  =  ^7  y(x)2  +(1  —  7r(x))2 

,  (l-7T(x))X-7  (x)y 

X0  (x)  = - r - r - 7- 

7  (x)“  +(1-7a.(x))" 


(20) 


(21) 


7v(x)2  +  (1  —  7a.(x))2 


//- 


/?oW'=- 


(1-7; r(x))  +  7v(x) 


•y'  '  F 


[7y(x)2  +  (l-7,(x))2 


(22) 


For  a  constant  velocity  or  uniform  motion: 


7*U)  =  lx  =  ~ 
v . 


(23) 


r  t  Vv 

7VU)=7V  =  — 


(24) 


and  m(x) ,  x0(x) ,  R0(x)  become  constant  m,  x:) ,  R()  in  (2)  defined  by  (20) -(22)  with 

7.r,yO)— >7x,y- 
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Using  (13),  the  system  Green’s  function  for  a  uniformly  moving  rigid  body  of  scatterers  can  be 
written  as: 


G( r, ux,X,Y,  jx,'Yy)  =  A( xQ  - 


RQuja 


r)' 


2km2[(l-°^(r-Rn^x)))2  ~  -^Y 

k  4km 

P(  —(r-Rn(cox)))exVU—(r-Rn(iox))2)- 

a  a 

/  2  I 

exp(-j2kRQ  [(1  -  ^-(r  -  Rn (wT)))2  -  ^  ■,  I2 ) exp(-jwvx0) 
k  4k  m 


(25) 


where  a'  ,Rq(lox)  are  defined  by  (ll)-(l 2)  and  constants  m,x0,R0  are  given  by  (20) -(22) 
(  7V  v  (x)  — >  7V  ).  Thus,  the  overall  range- Doppler  characteristic  of  a  moving  cluster  or 
distribution  is  described  by: 


’  D 


(uX’r’7x’7y)  =  ff  fD(X,Y,7x,7y)G(r,ux,X,Y,7x,7y)dXdY 


(26) 


where  fD(X,Y,  7^,7  )  represents  the  moving  scatterers’  stationary  scattering  distribution 
function.  Dependence  of  the  above  distribution  function  on  jx  and  jy  accounts  for  the 
variations  in  scattering  or  RCS  fading  due  to  the  movement.  If  these  effects  are  negligible,  the 
velocity  dependence  may  be  dropped,  i.e.,  fD(X  ,Y)  ,  and  the  motion  effects  will  be  included  in 
the  range -Doppler  Green’s  function. 
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3  Non-Uniform  Motion  Correction 


The  analytic  range -Doppler  frequency  model  developed  is  accurate  for  constant  velocity  target 
motion  assumptions.  The  method  can  be  generalized  to  non-uniform  or  varying  velocity  motions 
by  assuming  the  conservation  of  large  time-bandwidth  product  (TBP).  This  is  done  by 
introducing  a  method  for  correcting  stationary  point  of  the  variational  phase. 


The  mathematics  of  the  addressed  correction  technique  is  more  involved  due  to  the  nonlinear 
motion.  Here,  however,  certain  assumptions  are  made  that  are  judicious  and  can  provide  a 
feasible  yet  accurate  methodology.  The  stationary  point  equation  is  given  by  (5)  that  is  valid  for 
any  type  of  motion.  This  stationary  point  location  is  dependent  on  the  non-uniform  motion  due  to 
the  associated  range  function  variations,  i.e.,  left -hand-side  of  equation  (5).  As  addressed  in  the 
definition  of  (2),  for  non-uniform  motion  the  range  function  is  given  by: 


R(x)=  ^R()(x)2  +  m(x)2(x  -x0(x))2 


(27) 


where  m( x) ,  .^(x) ,  /^(x)  are  described  by  (20)-(22). 


The  exact  solution  of  (5)  for  the  corrected  stationary  point  can  be  written  in  a  general  form: 


Xs  =  xs(ux,u),) 


(28) 


As  pointed  out,  this  solution  is  involved.  However,  if  the  velocity  variations  are  not  severe,  one 
may  assume  m(x)  ,x0(x)  ,R0(x)  are  slow  varying  functions.  More  specifically  if: 


jWvf  d  d 

- - +  (x-x0(x))—  ln(m(x))-  —  ^(x)<l 

2m(x)(x—  x0(x))  dx  dx 


(29) 


upon  invoking  the  slow  varying  assumption  and  solving  for  the  stationary  points  (4)-(5),  one  can 
derive  the  corrected  Fourier  transform  as: 
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Sc (ux ,ur)  =  Cj  A(xq ( Xs (ux ,ur))~- 


_ RQ{Xs{ux,ur))ux 

2km2(X,(u>xMr)) +  , 


W, 


t)- 


2k  4k2m\X(ux,ur)) 


p(— ~~)exp(j  ~7~y 

2  a  4  a 


exp(-j2t«„(X,(u»J,u.,.)).[(l+^)2-  7 — - 

2k  4k  m  (X s(ux,ur)) 

exp  (~j'jjxx0(Xs  (ux,ur ))) 


i 

P). 


(30) 


To  derive  the  range-Doppler  transform  for  the  non-uniform  motion,  an  inverse  range  Fourier 
transform  on  (30)  is  to  be  applied  that  involves  the  stationary  point  solution  for  cur  ,  similar  to 
one  described  in  (7).  This  solution,  however,  is  more  complicated  due  to  dependence  of 
R{) ,  m  and  x0  functions  on  cor  .  In  order  to  develop  a  feasible  and  yet  accurate  solution,  a  two 
step  approach  is  undertaken  for  mathematical  manipulations.  As  the  first  step,  the  azimuth 
stationary  point  solution  (28)  to  the  general  stationary  point  equation  (5)  is  approximated  through 

(or  in  the  vicinity  of)  a  zeroth  order  solution  of  (7)  for  'jJr  .  In  the  next  step,  this  approximated 
stationary  point  is  used  to  define  the  overall  phase  function  in  (30)  to  solve  for  the  higher  order 
solution  of  (7)  for  ivr  .  Providing  that  the  exact  solution  of  the  stationary  phase  for  u)r  is 
variational  about  the  described  zeroth  order  solution,  the  associated  higher  order  solution  would 
be  accurate  reflecting  the  general  motion  due  to  the  variational  property  [11]-[12], 


To  obtain  the  zeroth  order  solution,  it  is  assumed  (further  to  the  slow-varying  assumption) 
functions  RQ(u>x,u)r)  ,m(u)x,u>r)  and  x0(ux,u)r)  are  variational  about  a  stationary  ur  point.  A 
judicious  choice  for  this  stationary  point  is  the  one  determined  by  the  average  motion  (i.e., 
velocity)  within  the  aperture  or  sub-aperture  of  interest.  Characteristics  of  this  point  are  defined 
by: 


m=  Vv  +(1_Tv)2 


(31) 


7/+(l-7*)2] 

h2  + 

X] 

(1-7*) +  7  yy 

2 

F2 

Tv2+a-Tv)2i 

(32) 
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(1-7JX-T yY 

7y2  +  a~%)2 


(33) 


where  and  jy  are  the  average  target  velocities  (relative)  within  the  aperture  of  interest. 
Accordingly,  similar  to  the  procedure  performed  in  (7)-(  10),  one  can  solve  for  the  zeroth  order 

•  •  —  f  cA 

stationary  point  ujr  using: 


2  a 


_ 

4  k2m2 


(34) 


where: 


a '  a 


.3-2  n  kV 


4>t  Vt7-(1  — 


4  k2  in2 


(35) 


Thus,  the  corrected  stationary  azimuth  point  (28)  can  be  approximated  by: 


xs  =  X  s(u>x,—2a'(r 


(1 


_A _ 

2  I 
_ \2 

4k2m2 


)) 


(36) 


where  Xs  (ux ,  u)r )  is  the  solution  of  (5)  for  the  general  target  motion. 


As  addressed  before,  the  higher  order  or  more  accurate  solution,  namely  jjJ v)  ,  for  the  stationary 

phase  can  be  found  by  using  (36)  to  define  the  overall  phase  function  in  (30).  Doing  so,  the  range 
frequency  stationary  point  satisfies  the  equation: 
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(37) 


^rzr+r-Ha(vx)=  0 

2  a 


where. 


a'  a 


Ro(xx)ux 


4k3m2(xs)(l- 


cu„ 


4k2m2(x  ) 


3 

)2 


Rfl(Ux)z 


R0(xs) 


(1- 


uv 


4k2m2(x  ) 


y- 


(38) 


(39) 


Applying  the  stationary  point  solution  (37)  to  perform  the  inverse  range  Fourier  transform,  one 
finds: 


Sc(u>x,r)=  CA(x0(xs)~ 


Ro(Xs)vx 


-12 


-)• 


P(  —(r-RQ(u;x)))exp(j—(r-Rn(ujx)Y)- 
a  a 

~r  2  i 

exp(— j'2fe/?0  (x  )[(\—°T-(r—  Rq  (ux  )))2  ] 2 )  exp(-  juxx0(xs )) 

k  4k  ”/?7 2 (x  ) 


(40) 


Equation  (40)  provides  a  relatively  rigorous  yet  analytically  feasible  means  to  model  non-uniform 
target  motion  in  the  range -Doppler  domain.  For  a  uniform  motion,  (40)  reduces  to  (13).  The 
stationary  point  (azimuth)  solution  in  (5)  appears  to  be  the  problematic  part  for  the  non-uniform 
motion.  In  most  cases,  however,  a  solution  can  be  found  through  analytic  derivations  and  proper 
approximations.  The  applied  approximations  yield  accurate  results  within  the  validity  of  the 
variational  and  slow-varying  assumptions.  These  assumptions  become  even  more  judicious  for 
sub-aperture  analysis.  For  piecewise  linear  (i.e.,  uniform  motion),  the  solution  approaches  the  one 
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in  (13)  within  the  associated  sub- aperture,  hence  providing  an  accurate  model  for  sub-aperture 
analysis. 
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4  Range-Doppler  Simulations 


The  effects  of  target  or  scatterer  motion  on  the  SAR  image  can  be  studied  by  the  aid  of  described 
range-Doppler  model.  In  doing  so,  motion  characteristics  may  be  investigated  through  the  whole 
aperture  or  at  the  sub- aperture  level.  The  latter  is  particularly  desired  when  characterizing  non- 
uniform  or  piecewise  uniform  target  motions.  In  practice,  the  measured  moving  point  target(s) 
range-Doppler  response  will  be  compared  against  the  simulated  model  data  to  extract  motion 
information  and/or  remove  the  associated  effects  from  the  SAR  image,  i.e.,  focusing.  Range- 
Doppler  characteristics  of  point  targets  with  uniform  and  non-uniform  motion  are  simulated  and 
presented.  Target  motion  parameters  and  SAR  imaging  system  specifications  are  assumed  to  be 
similar  to  those  of  a  typical  airborne  SAR  scenario  for  possible  further  validations.  Accordingly, 

it  is  assumed  that  the  flight  altitude  H  =  6  km  and  the  incident  angle  6inc  —  45°  . 


In  case  of  an  uniform  motion,  the  point  target  travels  at  a  constant  ground  speed  Vx  —  0. 1  va  and 
Vy  =  0.1  va  .  The  SAR  platform  or  vehicle  velocity  va  ,  is  typically  about  360  km/h  for  airborne 
SAR,  hence,  the  overall  ground  speed  of  the  moving  target  V-,,  =0.14  va  is  about  50  km/h.  For 
non-uniform  motion,  the  target  moves  with  a  similar  initial  speed  Vx  =  F  =0. 1  va  and 

accelerates  by  5  m/s2  at  each  direction.  Thus,  the  overall  ground  acceleration  reads  as 

9 

a  —  7. 1  m/s  •  For  non-uniform  motion,  the  analytic  range-Doppler  model  is  generalized  to 
accommodate  time  varying  velocity  by  correcting  the  point  of  stationary  phase.  Relatively 
accurate  Doppler  distribution  can  be  obtained  depending  on  the  validity  of  phase  variational 
assumption  that  is  good  for  typical  SAR  signals.  A  narrow  range  pulse  envelope  is  assumed  for 
current  simulations. 


Figure  3  depicts  the  non-uniform  motion  effects  on  the  range-azimuth  history  of  a  moving 
scatterer.  As  is  seen,  the  already  changed-shape  and  shifted  range  history  of  the  target  due  to  the 
uniform  motion  is  further  shape  changed  and  shifted  by  introducing  acceleration,  i.e.,  non- 
uniform  motion.  Thus,  the  target  image  will  be  further  smeared  and  shifted.  An  interesting 
observation,  however,  is  the  target  image  shift  that  is  partially  compensated  by  the  acceleration. 
Associated  range-Doppler  (instantaneous)  characteristics  of  the  described  movers  are  included  in 
Figure  4.  It  is  evident  the  non-uniform  mover  experiences  more  rapid  and  wider  Doppler 
variations.  Hence,  application  of  a  model  based  on  the  uniform  motion  assumption  within  the 
entire  aperture  can  not  be  effective.  Also,  it  is  seen  that  the  non-linear  motion  effects  are  more 
pronounced  at  the  higher  Doppler  frequencies.  Figure  5  illustrates  a  comprehensive  comparison 
amongst  the  range-Doppler  characteristics  derived  through  present  analytical  approach, 
instantaneous  using  the  exact  range  history,  and  instantaneous  using  the  Fresnel  approximated 
range  history  for  a  point  target  with  non-uniform  motion.  In  figures  4-5,  the  range-Doppler 
coordinate  system  (in  Doppler  direction)  is  centered  at  the  zero- Doppler  of  the  exact  range- 
Doppler  characteristic  of  the  uniformly  moving  target.  As  seen,  the  analytic  range-Doppler 
history  is  computed  applying  two  models.  One  is  based  on  the  exact  motion  model  for 
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generalized  velocity  functions  VJf)  and  Vy(t)  described  by  (15)-(16)  and  labeled  as 
Instantaneous.  The  other  is  based  on  the  same  equation  of  motion  but  with  velocities  defined  by 


the  instantaneous  derivatives  of  motion,  i.e.,  V  d  (t) 


dX(t) 

dt 


and  Vyder(t) 


dY(t) 

dt 


and  labeled 


as  derivative.  A  number  of  important  observations  are  made  through  the  inspection  of  Figure  5. 
Instantaneous  Doppler  variations  resemble  those  found  from  the  analytic  approach  using  the 
motion  derivatives.  This  is  somewhat  expected  as  instantaneous  Doppler  is  also  defined  by  the 
range  function  derivatives  that  involve  motion  derivatives.  Results  from  a  more  rigorous  analytic 
Doppler  characterization  (i.e.,  Inst.),  however,  differ  from  those  of  the  described  two.  The 
outlined  difference  (although  not  very  significant)  underscores  the  fact  that  the  instant  Doppler 
description  may  not  properly  and  accurately  represent  the  true  target  Doppler  distribution  for  non¬ 
linear  motion.  It  is  also  evident  that  Fresnel  phase  approximation  introduces  considerable 
difference  in  terms  of  range-Doppler  relations  for  moving  point  scatterer  compared  to  the  other 
techniques. 


Inspection  of  the  range  relation  (2)  and  system  Green’s  function  (25)-(26),  along  with  the  range- 
Doppler  simulations  of  point  scatterers  reveals  that  point  scatterers  belonging  to  a  moving  cluster 
or  target  do  not  theoretically  experience  the  same  Doppler  shift  that  results  in  shift  and  smearing 
in  the  image  plane.  Depending  on  the  severity  of  this  nonlinear  effect  that  is  related  to  a  number 
of  factors,  e.g.,  size,  resolution,  image  of  a  composite  target  may  be  further  distorted  and  difficult 
to  focus  using  conventional  methods  based  on  the  point  target  assumption.  If  so,  system  Green’s 
function  formalism  similar  to  (25) -(26)  can  be  applied. 
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Azimuth 

(r) 


Range 

(r) 


Figure  3.  Slant  Range -Azimuth  History  of  Point  Moving  Target  (yx  =  0.1, 7  =  0.1) 


cox  2k 


Range 

(r) 


Figure  4.  Instantaneous  Slant  Range-Doppler  History  of  Point  Moving  Target 

(7x  =  0.1,7  =0.1) 
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Analytic  range_Doppler(der.) 

-  Analytic  range_Doppler(inst.) 

.  Instantaneous  range_Doppler 

-  Fresnel_apprax  range_Doppler 

Figure  5.  Slant  Range -Doppler  History  of  Point  Moving 
Target(yx  —  0.1, yy  =  0.1, accel  =  a  =  7.1) :  Comparison 
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5  Velocity  Extraction 


One  characteristic  application  of  the  range -Doppler  model  of  moving  scatterers  is  the  velocity 
evaluation.  Such  velocity  information  and  motion  characteristics  are  used  for  advanced  detection, 
focusing  imagery  and  related  target  identification/recognition  applications. 


Technically,  one  can  apply  the  developed  model  or  Green’s  system  function  (25)  to  the  range- 
Doppler  history  of  a  uniformly  moving  scatterer  and  extract  the  velocity  information.  For  a  group 
or  cluster  of  moving  scatterers,  however,  one  (in  theory)  would  need  to  solve  the  system  integral 
equation  (26)  to  derive  the  motion  parameters.  Uniform  (i.e.,  constant  velocity)  motion  is  a  good 
assumption  for  a  variety  of  practical  applications  and  targets  (e.g.,  land  targets).  The  model  can 
also  be  applied  for  non-uniform  motion  at  the  sub-aperture  level  where  velocity  variations  are 
negligible  and  motion  is  piecewise  linear. 


Inspection  of  (25)  indicates  that  to  characterize  the  motion  of  a  persistent  point  scatterer  on  the 
target  of  interest,  the  associated  trajectory  in  the  range-Doppler  domain  must  be  tracked.  This 
trajectory  is  given  by: 


(41) 


Accordingly,  the  moving  scatterer  path  in  the  range-Doppler  plane  is  characterized  by: 


r 


(1 


_ 

2  I 

Ak2m2 


(42) 


Using  (20)-(24),  relation  (42)  may  be  written  as: 


(m2H2+Y2B2)2 

r  =  - - (43) 

where 

m2  =  7y  2  +  (1  —  7_r  )2  (44) 
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B  =  (  l-7,)+7rf 


(45) 


and  £  =  — is  the  normalized  Doppler  frequency. 
2& 


According  to  the  above,  the  problem  of  moving  scatterer  tracking  and  velocity  extraction  would 
effectively  transform  to  the  optimal  velocity  (~/x,  7  )  derivation  for  the  respective  range-Doppler 

history  function.  A  least -square  based  technique  is  applied  to  compute  the  optimal  velocity 
quantities  from  the  range-Doppler  history.  However,  instead  of  performing  the  optimization  in 
terms  of  the  relative  velocity  parameters  (jx,  7 v) ,  the  error  target  function  is  minimized  with 

respect  to  functions  (44) -(45): 


(46) 

B=fl(  7,.7y) 

(47) 

This  can  be  shown  to  provide  a  more  effective  optimization. 

Define  the  target  function: 

(/i2(7„7v)^2+^2/22(7„7y))2 

y=F{fl{ix,iy)J2iix,iy),0=K- - - ! — — 

(48) 

[fAlx’ly)  -£2)2 


Accordingly,  the  least-square  optimization  scheme  is  governed  by  the  equation: 


4  =  1,2  (49) 


with  "i"  and  "  N"  denoting  the  sample  index  number  and  total  number  of  samples  within  the 
aperture  or  subaperture  of  interest. 
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As  addressed,  the  choice  of  /j  (7^,7  ),f2  (7*>7V)  yields  analytic  solution  to  (49)  for  the  relative 
velocities  77,77  .  Analytic  derivations  and  solutions  are  included  in  Annex  A. 


In  order  to  evaluate  the  efficiency  of  the  model  and  the  velocity  extraction  algorithm,  one  should 
apply  the  algorithm  to  the  range-Doppler  history  of  a  moving  scatterer  with  known  motion.  To  do 
this,  the  range-Doppler  characteristics  of  a  point  scatterer  moving  with  constant  velocity  are 
simulated  and  used  for  velocity  extraction  applying  the  described  procedure.  These  derivations 
and  velocity  estimations  are  performed  for  two  cases  of  uniform  motion,  i.e.,  jx  —  7  =0.1  and 

')x  =  0. 17 ,7  =  0.23  .  These  motion  parameters  are  chosen  to  somewhat  emulate  two  typical  and 
diverse  cases  of  motion  and  operation. 


The  range-Doppler  history  simulations  for  the  above  two  cases  in  the  range-azimuth  and  range- 
Doppler  frequency  domains  are  depicted  in  Figures  6-9.  The  motion  parameters  used  in 
simulations  and  respective  computed  velocities  are  included  in  Tables  1,  2.  As  seen,  a  higher 

incident  angle  dinc  =53.13°  is  chosen  that  is  more  typical  for  scenarios  including  ground  moving 
targets.  Moreover,  to  make  it  more  general,  a  non- broadside  moving  target  (i.e.,  non- zero  Doppler 
at  the  initial  state  t  =  0  )  is  considered.  The  slant  range-azimuth  history  of  this  non-broadside  point 
target  is  presented  in  Figure  6  for  a  better  understanding  of  the  non-zero  Doppler  (i.e.,  squint) 
imaging  of  moving  targets. 
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Range 

(r) 


9.9698E3  9.9699E3  9.97E3  9.9701E3  9.9702E3  9.9703E3  9.9704E3 


Figure  6.  Slant  Range-Azimuth  History  of  Point  Moving  Target  (yx  =0.1, 7,,  =  0.1,  X  =100  in) 


Figure  7.  Instantaneous  Slant  Range-Doppler  History  of  Point  Moving  Target 
(7a  =0-l>7y  =  0.1,X  =100m) 
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Table  1  includes  computed  results  for  velocity  extraction  using  the  simulated  range-Doppler 
history  of  a  mover  with  motion  parameters  7^  =  7  =0.1  (Figures  6-7)  in  10  different 

subapertures  within  the  total  aperture.  The  sizes  of  subapertures  (Doppler  width)  are  chosen  to  be 
equal  while  utilizing  the  whole  available  aperture,  i.e.,  the  available  aperture  divided  into  10 
equal  segments.  As  can  be  seen,  computed  velocity  parameters  are  in  good  agreement  (less  than 
1  %  error)  with  the  parameters  used  in  simulations  or  equivalent  truthing. 


Motion  Parameters  Used  in  the  Range 
History  Simulation 

Extracted/Calculated  Motion  Parameters 

Subaperture 

# 

lx 

ly 

H(m) 

N 

2 

m 

B 2 

lx 

ly 

1 

0.1 

0.1 

100 

8000 

6000 

10 

.820415 

.812662 

.099772 

.100024 

2 

0.1 

0.1 

100 

8000 

6000 

10 

.820198 

.812448 

.099891 

.100011 

3 

0.1 

0.1 

100 

8000 

6000 

10 

.819896 

.812149 

.100056 

.099993 

4 

0.1 

0.1 

100 

8000 

6000 

10 

.819494 

.811750 

.100277 

.099969 

5 

0.1 

0.1 

100 

8000 

6000 

10 

.819905 

.812158 

.100051 

.099994 

6 

0.1 

0.1 

100 

8000 

6000 

10 

.820510 

.812757 

.099719 

.100031 

7 

0.1 

0.1 

100 

8000 

6000 

10 

.819724 

.811978 

.100151 

.099983 

8 

0.1 

0.1 

100 

8000 

6000 

10 

.820883 

.813126 

.099515 

.100053 

9 

0.1 

0.1 

100 

8000 

6000 

10 

.820090 

.812340 

.099950 

.100005 

10 

0.1 

0.1 

100 

8000 

6000 

10 

.820122 

.812373 

.099932 

.100007 

Table  1.  Simulated  and  Computed  Motion  Parameters  for  Point  Moving  Target 
(7a.  =  0.1, 7^  =0.1,X  =100/n) 
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Figure  8.  Slant  Range-Azimuth  History  of  Point  Moving  Target 
(7a-  =  0.17,  jy  =0.23,  X  =100 m) 


Figure  9.  Instantaneous  Slant  Range-Doppler  History  of  Point  Moving 
Target  (,-fx  =  0.17,7^  =  0.23,  X  =  100m) 
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Computed  velocity  parameters  for  a  mover  with  motion  parameters  jx  =0.17,7  =0.23  are 

provided  in  Table  2.  Simulated  range- azimuth/range-Doppler  characteristics  are  given  by  Figures 
8-9.  Again,  a  good  agreement  with  the  mover’s  true  motion  parameters  (simulated)  is  evident 
(less  than  1%  error).  To  examine  the  applicability  of  the  present  algorithm  to  various  parts  of  the 
available  Doppler  spectrum  (synthetic  aperture)  or  subapertures,  an  additional  subaperture  #11  is 
considered  for  velocity  extraction.  This  subaperture  is  arbitrarily  chosen  within  the  total  aperture 
and  has  a  different  width  size  and  center  location  from  other  subapertures. 


Motion  Parameters  Used  in  the  Range  History 
Simulation 

Extracted/Calculated  Motion  Parameters 

Subaperture 

# 

Tx 

ly 

X(m) 

Y(m) 

H  (m) 

N 

2 

m 

B2 

lx 

ly 

1 

0.17 

0.23 

100 

8000 

6000 

10 

.741641 

.693531 

.170089 

.229975 

2 

0.17 

0.23 

100 

8000 

6000 

10 

.741784 

.693666 

.170008 

.229997 

3 

0.17 

0.23 

100 

8000 

6000 

10 

.741789 

.693671 

.170005 

.229998 

4 

0.17 

0.23 

100 

8000 

6000 

10 

.741682 

.693570 

.170065 

.229981 

5 

0.17 

0.23 

100 

8000 

6000 

10 

.742027 

.693893 

.169872 

.230035 

6 

0.17 

0.23 

100 

8000 

6000 

10 

.742027 

.693893 

.169872 

.230035 

7 

0.17 

0.23 

100 

8000 

6000 

10 

.741606 

.693499 

.170108 

.229969 

8 

0.17 

0.23 

100 

8000 

6000 

10 

.741514 

.693414 

.170159 

.229955 

9 

0.17 

0.23 

100 

8000 

6000 

10 

.741784 

.693666 

.170008 

.229997 

10 

0.17 

0.23 

100 

8000 

6000 

10 

.741725 

.693611 

.170041 

.229988 

11 

0.17 

0.23 

100 

8000 

6000 

20 

.742177 

.694033 

.169789 

.230058 

Table  2.  Simulated  and  Computed  Motion  Parameters  for  Point  Moving  Target 
(7,  =  0. 17, 7y  =0.23,  X  =  100m) 
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6  Conclusions 


The  appropriate  choice  of  the  backscattering  signal  model  is  highly  important  for  effective  SAR 
data/image  exploitations.  The  need  for  such  a  robust  analytic  modeling  is  more  critical  for 
operations  involving  moving  targets.  There  exist  conventional  models  with  their  own  set  of 
advantages  and  shortcomings  depending  on  the  type  of  exploitation  applications.  Considering  the 
coherent  principle  and  nature  of  SAR  imaging  and  associated  products,  accuracy  of  the  SAR 
signal  model  is  provided  through  accurate  characterization  of  the  basic  spherical  phase.  On  this 
basis,  an  analytic  range-Doppler  model  to  characterize  canonical  moving  point  scatterers  is 
developed.  This  coherent  model  is  used  to  establish  the  system  response  or  Green’s  function  for 
imaging  a  collection  or  distribution  of  moving  scatterers.  The  approach  is  accurate  as  the 
spherical  phase  front  variations  and  Doppler  effects  due  to  the  target  motion  are  directly  modeled 
in  slant  range-Doppler  domain  without  the  need  for  range  and  azimuth  de -coupling.  The  slant 
range-azimuth  Green’s  function  formalism  is  mathematically  robust  as  it  utilizes  the  signal  space 
data,  and  reflects  the  electromagnetic  physics  of  scattering  as  the  exact  Green’s  function  for  a 
point  source  is  defined  in  slant  range.  Also,  the  range  band-limited  property  is  incorporated  in  the 
formulation  of  target  Doppler  distribution. 


The  developed  model  enjoys  certain  attractive  characteristics.  Accurate  modeling  of  range 
variations  and  its  azimuth  dependence  provides  better  prediction  of  range-Doppler  behaviour  of 
the  moving  scatterer.  Higher  squint  angle  and  wider  synthetic  aperture  modes  of  operation  would 
benefit  from  such  accuracy,  where,  conventional  phase  expansions  (e.g.,  series,  Fresnel)  would 
not  be  suitable.  Furthermore,  the  target  dependent  shift- varying  property  (conventionally  ignored 
through  the  plane  wave  approximation  process)  is  captured.  Application  of  the  range  -Doppler 
model  in  conjunction  with  the  measured  data  or  image  yields  target  velocity  information.  It  can 
also  be  used  to  define  the  matched-filter  transfer  function  for  image  re-focusing.  The  means  for 
sub-aperture  analysis  is  provided  which  enables  characterization  of  motion  variations  within  the 
synthetic  aperture  acquisition  time.  Results  shown  indicate  accurate  derivation  of  velocity 
parameters  for  a  given  mover  using  various  subapertures  of  range-Doppler  data  (in  terms  of 
center  location  and  width  size).  The  model  is  generalized  to  treat  nonlinear  target  motion.  This 
generalized  model  is  relatively  involved  but  is  accurate  and  practical  to  study  and  treat  the 
localized  effects.  Moreover,  adopting  a  sub- aperture  approach  and  assuming  piecewise  linear 
motion,  the  uniform  motion  range-Doppler  model  can  be  locally  applied. 


The  key  contribution  of  the  present  work  may  be  viewed  as  providing  a  canonical  picture  of 
motion  effects  on  the  range-Doppler  frequency  characteristics.  One  important  observation  is  that 
the  scatterer  motion  is  more  consistently  analyzed  and  modeled  in  slant  range-Doppler  domain. 
Accordingly,  the  model  validation  and  required  experimentation  are  more  meaningful  in  such 
domain. 


Actual  experimentation  based  on  real  imagery  will  be  presented  in  an  ensuing  communication  due 
to  unavailability  of  such  results  at  the  time  of  submission.  Furthermore,  considering  the 
importance  of  experimental  verifications  a  complete  paper  ought  to  be  dedicated  to  the 
experimentations,  validations  and  interpretations.  There  are  no  standard  measures/metrics  for 
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focusing  evaluation  that  are  solely  based  on  the  final  focused  image.  More  objective  assessments 
incorporating  the  focused  target  response  signal  characteristics  or  at  the  intermediate  stage  (e.g., 
range -Doppler)  is  needed  for  proper  validations.  Such  measured  validations  become  more 
important  where  the  imagery  products  are  used  for  classification  and  recognition  purposes. 
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Annex  A  Velocity  Analytic  Description 


Solution  to  the  range -Doppler  error  target  function  optimization  (49)  for  moving  scatterer 
velocity  components  can  be  derived  through  considerable  amount  of  analytic  manipulations  (to 

'3  £4  ' 


3 

the  order  of  - -A- ,i.e..O 
8  fi 


8  m 


),  as: 


m^=fl2(7x,7y)=i- 


,  E*E$4-EtfE*$2 

A  AT  AT  AT  AT 


2  *E*i -E*E$a 


(Al) 


B2=f22( 7x,7y)= 


a  /r  2  2  tt  2 

M  —m  H 


(A2) 


where 


w5E^+Am3E^2 

_ N _ Z  N _ 

Nm4  +  m2J2^+]AJ2^ 

N  4  N 


(A3) 


Velocity  components  7^ ,  7y  can  be  found  by  solving  the  system  of  equations  (44)-(45)  using  the 
results  for  m,B  obtained  from  (A1)-(A3). 
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DND 

Department  of  National  Defence 
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Defence  Research  &  Development  Canada 

DRDKIM 

Director  Research  and  Development  Knowledge  and  Information 
Management 

ISAR 

Inverse  Synthetic  Aperture  Radar 

PSP 

Principle  of  Stationary  Phase 

PolSAR 

Polarimetric  Synthetic  Aperture  Radar 

R&D 

Research  &  Development 

SAR 

Synthetic  Aperture  Radar 

TBP 

Time  Bandwidth  Product 
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